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We examine the dynamics of the gravitational collapse in a 4-dim Lorentzian brane em- 
bedded in a 5-dim bulk with an extra timelike dimension. By considering the collapse 
of pure dust on the brane we derive a bouncing FLRW interior solution and match it 
with a corrected Schwarzschild exterior geometry. In the physical domain considered for 
the parameters of the solution, the analytical extension is built, exhibiting an exterior 
event horizon and a Cauchy horizon, analogous to the Reissner-Nordstrom solution. For 
such an exterior geometry we examine the effects of the bulk-brane corrections in the 
Hawking radiation. In this scenario the model extends Bekenstein's black hole geomet- 
rical thermodynamics for quasi-extremal configurations, with an extra work term in the 
laws associated with variations of the brane tension. We also propose a simple statisti- 
cal mechanics model for the entropy of the bouncing collapsed matter by quantizing its 
fluctuations and constructing the associated partition function. This entropy differs from 
the geometrical entropy by an additive constant proportional to the area of the extremal 
black hole and satisfies an analogous first law of thermodynamics. A possible connection 
between both entropies is discussed. 
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1. Introduction 

Black holes are solutions of vacuum general relativity equations describing the ex- 
terior spacetime of the final stage of gravitationally bounded systems whose masses 
exceeded the limits for a finite equilibrium configuration!^ Geometrically a black 
hole may be described as a region of asymptotically flat spacetimes bounded by an 
event horizon hidin g a singularity formed in the collapse. Fundamental theorems 
by Israel and Cartei^El state that the final stage of a general collapse of uncharged 
matter is typically a Kerr black hole, which has an involved singularity structure. 
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Nevertheless, for a realistic gravitational collapse we have no evidence that the 
Kerr solution describes accurately the interior geometry of the black hole. On the 
contrary, the best theoretical evidence presently available indicates that the interior 
of the black hole thus formed is analogous to the interior of a Schwarzschild black 
hole with a global spacelike singularity.^ The simplest way of forming such structure 
is by the spherical collapse of dust, as originally shown in the classical paper of Op- 
penheimer and Snyder!^ However, as singularities cannot be empirically conceived, 
this turns out to be a huge pathology of the theory. 

Notwithstanding the cosmic censorship hypothesis^ (CCH), there is no doubt 
that the general theory of relativity must be properly corrected or even replaced 
by a completely new theory, let us say a quantum theory of gravity. This demand 
is in order to solve the issue of the presence of singularities predicted by classi- 
cal general relativity, either in the formation of a black hole or in the beginning 
of the universe. While a full quantum gravity theory remains presently an elusive 
theoretical problem, quantum gravity corrections near singularities formed by grav- 
itational collapse have been the o bje ct of much recent research, from loop quantum 
cosmology-^ to D-brane theory!Srll2J j n the latter scenario extra dimensions are in- 
troduced constituting the bulk space. All matter would be trapped on a 4-dim 
world-brane spacetime embedded in the bulk and only gravitons would be allowed 
to move in the full bulk. At low energies general relativity is recovere d 8 but at high 
energy scales significant changes are introduced into the gravitational dynamics and 
the singularity could be eventually removed. 

The problem of the gravitational collapse in the bra neworld scenario has been 
the object of several important works. Bruni et al\^ studied the Oppenheimer 
Snyder collapse on a Randall-Sundrum-type brane and showed that the exterior 
vacuum spacetime on the brane cannot be static and there fore precluding the for- 
mation of black holes in the theory. However Dadhich et aZffSl have demonstrated 
the existence of static black holes on the brane in the Randall-Sundrum scenario. 
These black holes are exact solutions of the effective Einstein equations on the brane 
and correspond to ReissnerNordstrom (RN)-type black holes, with a tidal charge 
(also denoted KaluzaKlein (KK) charge) originated from the 5-dim Weyl curvature 
instead of an electric charge. In this vein, Govender and DadhiclP^ constructed a 
model of the OppenheimerSnyder collapse in the brane in which the collapsing solu- 
tion is matched to the brane generalized Vaydia solution which in turn is matched 
to the asymptotically flat RN-type metric with a KK charge. The mediation by 
the Vaydia radiation metric is a new feature introduced by the Randall-Sundrum- 
type brane so that the collapsing sphere radiates null radiation. This picture is 
the paradigm of the gravitational collapse of a homogeneous spherically symmet- 
ric configuration on a Randall-Sundrum-type brane embedded in a nonconformally 
flat, but otherwise vacuum bulk. The problem of the gravitational collapse of a null 
fluid on the brane was approached by Dadhich and Ghoshp^- 1 where the parameter 
windows in the initial data set, giving rise to a naked singularity or favoring the 
formation of black holes, are examined. 



December 12, 2012 2:7 WSPC/INSTRUCTION FILE text-ijmpd 



Gravitational Collapse and Black Hole Thermodynamics in Braneworld Scenario 3 

In our approach in this paper we have considered the gravitational collapse of 
a spherically symmetric dust distribution in the framework of a braneworld sce- 
nario, with a 5D bulk having an extra timelike dimension. The interior geometry 
is still given by the FLRW metric but - due to the timelike character of the extra 
dimension - the dynamics of the collapsing dust has an effective potential barrier 
generated by the bulk-brane corrections. This potential barrier avoids the formation 
of a singularity yielding a perpetually oscillating collapsed matter. We obtain the 
unique static exterior geometry which is smoothly matched to the interior FLRW 
geometry. 

As we know from General Relativity, the CCH addresses the issue whether a sin- 
gularity thus formed in gravitational collapse is visible to an asymptotic observer 
or hidden by an event horizon. As we shall see in our model, if the total mass of 
the collapsing dust is larger than a critical value, we obtain a static solution which 
corresponds to a nonsingular black hole with an event horizon (besides a Cauchy 
horizon) encapsulating not a singularity but a perpetually oscillating collapsed mat- 
ter. Although this perpetually bouncing matter is not visible from an asymptotic 
observer, no singularity is engendered. This new feature gives rise to a modified 
CCH which now addresses the issue whether the perpetually bouncing matter is 
visible to asymptotic observer or hidden by an event horizon. 

By considering our exterior static solution, we construct a statistical model for 
the quantum degrees of freedom of the oscillating collapsed matter whose entropy 
can be associated with the entropy fluctuations about the extremal configuration of 
the exterior static black hole geometry. In this direction we are also led to evaluate 
the Hawking evaporation processes of the exterior black hole. The no-go theorem 
of Ref. 14 may be circumvented as long as one relaxes the condition of a vacuum 
nonconformally flat bulk assuming, for instance, a bulk matter content sa tisfy ing 
energy conditions or the presence of torsion degrees of freedom in the bulkU^ We 
have however not addressed the 5D equations for the determination of the bulk 
space since this task is beyond the scope of this paper. 

For the sake of completeness let us give a brief introduction to braneworld theory, 
making explicit the specific assumptions used in obtaining the dynamics of the 
model. We rely on Refs. 9-12, and our notation basically follows Let us start 
with a 4-dim Lorentzian brane S with metric g a b, embedded in a 5-dim conformally 
flat bulk M with metric gAB- Capital Latin indices range from to 4, small Latin 
indices range from to 3. We regard S as a common boundary of two pieces Mi 
and M.2 of M and the metric g a b induced on the brane by the metric of the two 
pieces should coincide although the extrinsic curvatures of E in jVfi and M2 are 
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allowed to be different. The action for the theory has the general form 



1 

24 



r(5) 



Mi 



R - 2A 5 

(5) 



dfx 



R - 2A 5 



2n\L 5 



d 5 x + 2e / y/-^)gK 2 d' l x 



-2c J sT^gK^x) + i jf ^Wg( K ^ ) R-2a)d i x 



+ / V^~gL 4 (g a ^p)d 4 x. (1) 



In the above ^R is the Ricci scalar of the Lorentzian 5-dim metric gAB on M, 

and ^ R is the scalar curvature of the induced metric g a b on E. The parameter a is 

denoted the brane tension. The unit vector n A normal to the boundary E has norm 

e. If e = —1 the signature of the bulk space is (—,—,+,+,+), so that the extra 

dimension is timelike. The quantity K = K a b g ab is the trace of the symmetric 

tensor of extrinsic curvature K„h = Y^ a c Yj, D ^c n D, where Y A (x a ) are the 

embedding functions of E in While L^g^.p) is the Lagrangean density of 

Jonl 

the perfect fluid- - (with equation of state p = ap), whose dynamics is restricted to 
the brane E, L5 denotes the lagrangian of matter in the bulk. All integrations over 
the bulk and the brane are taken with the natural volume elements y— e^hj d 5 x 
and y — ( 4 )g d 4 x respectively. K5 and K4 are Einstein constants in five and four- 
dimensions. With the exception of Section 6, throughout the paper we use units 
such that K = c = 1. 

Variations that leave the induced metric on E intact result in the equations 



(5), 



>G AB +A 5 Wg AB =K^ 5 >T AB , (2) 
while considering arbitrary variations of g AB and taking into account ([2]) we obtain 



(4) G Qb 



«5 



\ D ab 



s 



K4 



[jab — crg a bj , 



(3) 



where S a b = Kab — Kg a b- In the limit K4 — > 00 equation © reduces to the Israel- 
Darmois junction conditional] 



S 



(i) 



(jab 



og a b 



(4) 



We impose the ^-symmetrjH^ and use the junction conditions (U]) to determine 
the extrinsic curvature on the brane, 



K ab = 



:«5 



(T a b 



1 , O 

2 T 9ab) + g gab 



Now using Gauss equation 

{4) Rabcd = {5) R M NB.sY^Y^Y^Y s d + e(K ac K hd - K ad K bc ) 



(5) 



(6) 
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together with equations @ and (|S|) we arrive at the induced field equations on the 
brane 

^G ab = -A 4 (4) 3ab + 87rGATT ab + e4n ab - eE ab + eF ab (7) 
where we define 

A 4 :=^(| + ^ 2 ), (8) 

G N :=e±, (9) 
Ha, := ~\r c a r bc + ±TT ab + l^g ab r cd r cd - ^r^g ab , (10) 
F ab := ^T BD Y B a yf + [( 5 )T BI3 n B n D - i e ^t] ( 4 ). 9ah }. (11) 

Ga? is just the Newton's constant on the brane. Here we remark that the effective 
4-dim cosmological constant can be set zero in the present case of an extra timelike 
dimension, by properly fixing the bulk cosmological constant as A5 = i«| a 2 . It's 
important to notice that for a 4-dim brane embedded in a conformally flat bulk we 
have the absence of the conformal tensor projection and F ab in Eq. (|8|). Accordingly 
Codazzi's equations imply that 

V a K-V b K b a = ~e K 2 5 V b r b a , (12) 

By imposing that V(,r„ = the Codazzi conditions read 

V a E ab = 4V a U ab + V a F ab . (13) 

where V a is the covariant derivative with respect to the induced metric g ab . Equa- 
tions (0 and ([T3J) are the dynamical equations of the gravitational field on the 
brane. 



2. The Interior Solution and the Exterior Geometry 

We assume a spacetime braneworld model embedded in a 5-dim dc Sitter bulk with a 
timelike extra dimension (e = —1), whose matter content is a spherically symmetric 
collapsing dust with density p. In a coordinate system comoving with d ust the 
interior geometry is still shown to be a Friedmann-Robertson- Walker metriJ^ 

ds 2 = -dt 2 +a 2 (t)( l — dr 2 + r 2 dn 2 ) (14) 

V 1 — kr z ) 

with its dynamics given by the first order modified Friedmann equation 
.3 , , 8nG N E 4nG N E 2 

a =-«H 5 Q | 1 4 (15) 

where we made A4 = by a proper choice of A5. Eq is a constant of motion 
associated with the dust density, Eq — pa 3 , and a is the negative brane tension so 
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that Gn > 0. Assuming initial conditions for the collapse d(0) = and a(0) = 1, 
we get 

3 L ° 2 a 



(16) 



By defining 



Pa 



a , V(a) := 



2ttG n E$ 4ttG n E ( 



N&0 



3\a\a 4 

one may infer the Hamiltonian constraint 



3o 



H 



Pi 



+ V(a) = 0, 



(17) 



(18) 



which is equivalent to the first integral (|15|) . An immediate calculation shows that 
from this constraint one can easily obtain the motion equation 

aa + 2a + 2k = . 

a 

From now on we are going to assume k > and that the potential V(a) has two 
real positive roots. These assumptions restrict the domain of the parameters as 



\a\ > 2E Q . 



(19) 



Therefore we see that V(a) has one extremal located at a = (^j^fj > and two 
positive real roots a\ — a m i n and ei2 = 1 (with a m .;„ < 1). According to the 
restriction (TT9|) . the potential V(a) gives us an oscillatory solution for the scale 
factor (between a m i„ < a < 1) avoiding the singularity formation at the center of 
the matter distribution at a = 0. 

Let us now consider the following coordinates transformation 



R = ar , 6 = 6, (f> = (f>. 
In this sense, the line element can be written as 



ds z = - 



a 2 R 2 



1 - 



kR 2 



dt 1 - 



2aaR 
i 2 - kR 2 



dtdR 



a 2 - kR 2 



(20) 



dR 2 + R 2 dn 2 , (21) 



where dft 2 = dil 2 . Defining 

T := F[S(R, t)], C(R,t):= 



1 - 



a 2 R 2 



we get: 



E(R, t) 



ds =- c (m) (is) dT +2 (m) (ds) Het) (or) +e 



aaR 



i 2 - kR 2 

„2 



kR 2 



D(R,t) 



'dS\- 1 sdF\-i 



kR 2 ' 



'dS\- 1 sdS 



dt) KdRJ 



dt) \dRJ 



dTdR 
dR 2 + R 2 dfl 2 . (22) 
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According to the Birkhoff theorem^] j n General Relativity, we know that 
the exterior solution of a spherically symmetric collapse of dust is given by the 
Schwarzschild geometry where the metric is diagonal. As we are motivated to find 
a correction of the Schwarzschild geometry, let us consider the following condition 



C 



\dt) \dRJ 



E = 0. 



Therefore the metric (|22|) is diagonal and we automatically guarantee that 



9RR = 



td2 

l_iL(fc + a 2 



It is easy to verify that the solution for (|2"3")l is given by 

S(R, t) = S+ ^^a 2 -kR 2 exp ( - k [ —, 
a \ J aa< 

where 5 and \i are arbitrary constants, 



da 



(23) 
(24) 
(25) 



;da = Din \a — 1| — /In \a - 
2J 



with 



A. . 1 -\- &min 

D := 
I :-- 



^AB - A 2 



arctan 



(■ 



3k 



B := 



2a + A 
V \/AB - A 2 ' 

^G N E 2 
3k\<7\o Jm i n 



(26) 



1 



{l-a min )(A + B + 1) 
a 3 

min 

(1 — 0>min)(ami7i + 

mm ~r J-J j 



B(Da r , 



E) 



J 

Let us now assume that the junction of the interior solution with the exterior 
geometry is given at the surface defined by r = 7, where 7 determines the boundary 
of the matter distribution. By defining the constant M := 47r7 3 E ° , we obtain 



9RR 



1 



2G N M 3G N M 2 
_ 



(27) 



Air\a\R 4 . 

Employing the integrating factor technique J22l we define the function F[S(R, t)] 
by the following differential equation 

dF_ 
~dS 



,2* r^T 



kR 2 



(28) 



i 2 -R 2 {k + a 2 ) L k(S-S) 

where t and R are implicit functions of S through equation (|25[) . Adopting this 
choice we verify from equation (|22[) that 

1 



grr 



R 2 

1 - ^(fc + a 2 ) 



9RR 



(29) 
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It is important to remark that when r = 7, S(a) is a monotonous function in 
the physical domain of aQ which may be properly inverted in such a way that we 
can express a in terms of S. 

Therefore, from equations (|27|) and (|29|) we see that the exterior geometry reads 



ds 2 = - 



1 - 



2G N M 3G N M 



R 



4 7 r|cr| J R 4 



dT 2 



1 - 



2G N M 3G N M 



R 



4 7 r|er| J R 4 
-R 2 (d0 2 + 



dR 2 



'■)■ (30) 



From (|15p and (|30l) we can see that General Relativiy is recov ered, with the 
corresponding standard Oppenheimer-Snyder models, when \a\ — > 

In the following Sections we discuss some properties of the exterior geometry 
()30[) and construct its maximal analytical extension in order to better understand 
the avoidance of the Schwarzschild singularity. We also discuss some physical issues 
in this geometry as the Hawking temperature of the black hole and an extension of 
black hole thermodynamics in this spacetime. 



3. Analytic Completion of the Manifold 

In order to examine the analytic completion of the exterior geometry (|30j) , we need 
to know whether, and under what circumstances, the configuration forms event 
horizons. By defining the polynomial P(R) ■— R^grr, we see that a necessary 
condition for horizon formation is that the mass M, of the collapsing star, equals 
or exceeds a critical limit M* , 

M ^ M *=(9^r< (3i) 

(\a\ fixed), namely, that the polynomial P(R) has one i?*, or two < R + ) roots 
respectively (cf. Fig. 1). Otherwise we cannot have formation of event horizons. For 
illustration let us consider again the parameters given in (a) but now with a large 
value of the brane tension (|ct| = 10 8 ) to tentatively approach general relativity. It 
is easy to check that we obtain again a strictly monotonous function S(a) in the 
physical range of a, and i?_ > ^a m in ~ specifically 7a m j„ = 0.02597827765 and 
i?_ = 0.02597828531. As \a\ decreases, the difference —ja m i n ) increases by the 
same order of magnitude. 

On the other hand, we note that the critical mass depends solely on the parame- 
ter \a\. To have an idea of the order of magnitude of | cr | for event horizon formation, 
let us take M* ~ 1.4M©, the Chandrasekhar limit. This yields |cr| ~ 10~ 2 fcm -2 . 
A star with the Chandrasekhar mass will not form an event horizon if the brane 



a In order to illustrate the behavior of the integral H26H in the physical domain of our parameters, 
let us assume |ct| = 0.05m -2 , Gjv = 1, E = 1.04 x lO" 23 ™ -2 , 7 = 6.96 x 10 8 m, 5 = 1, fl = 0.01. 
We are assuming a reasonable physical value for the brane tension while considering a density one 
hundred times greater than the density of the sun. By taking these parameters we see that S(a) is 
a monotonous function of a in the physical domain a m i n = 4.702669449 X 10~ 8 < a < a max = 1. 
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tension is smaller than 10 2 km 2 , so that this value establishes a lower bound for 

Considering then the case of two positive real roots, the polynomial P(R) may 
be rewritten as 



P(R) = (R- R+)(R - R-)(R 2 +aR + (3), 



where 



a = R+ + R_ - 2G N M , f3 = 



3G N M 2 



(32) 



(33) 



ATr\a\R + R- 

In this case the collapse of the surface of dust must cross i?_ , that is, ja m < R- so 




Fig. 1. The polynomial P(R) = R 4 gTT for dust masses M < M„ (no event horizon), M = M* 
(extremal black hole) and M > M* (black hole formation with a event horizon and a Cauchy 
horizon R-). We consider here \a\ = 0.05, in units Gjv = c= 1. 



that a stable black hole forms with trapped perpetually bouncing matter. 
Let us then consider the following coordinates transformation 



u R\R) 



2k , R 4 
— av := — ; — - 

v P{R) 



dR + dT. 



Therefore, from we get 



, 2 4k 2 P(R) 
ds — r—dudv 

uv R 4 



R 2 {d9 2 



(34) 
(35) 



December 12, 2012 2:7 WSPC/INSTRUCTION FILE text-ijmpd 



10 R. Maier and I. D. Soares 

By defining 




Fig. 2. Kruskal diagrams showing the curve ABCD which illustrates how the matter surface 
evolves in time. The shaded portion corresponds to the interior of the matter distribution. Figures 
(a) and (b) overlap in the region //(,, and may be regarded as linked together along the curve 
R = Rq , where Rg is any convenient value between R+ and R _ . 



R* = R + G\n 



where 





+ Wki 




JL_i 


±-1 


R+ 




R- 



+ 



2 J 



; arctan 



2R + a 



f- 



G 
W 
J 



R\ 



(R+ - R-)(R 2 + + aR+ + p) 



Ri 



(R+ -R-)(R 2 _ 

r>{ A B 

p(r- + + r-- 1 



Therefore, integration of (f34|) and (|35|) yields 

R* = k\ti\uv\ , T = k\tl\v/u\. 



(37) 

(38) 
(39) 
(40) 

(41) 



Consider now the chart (u±, v-y) (defined for R > obtained by setting n = G. 
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From d33) and ^ we find that 

/ R \{R 



x exp 



R 2 J / 2R 

arctan 



G Gy/4/3 - a 2 K VW 



and (|36p exhibits no singularity at R = The chart (ui, v\) in fact gives a regular 
mapping of any given subregion of the manifold which has R > R-. However, a 
coordinate singularity does develop at R = i?_ and it is necessary to go over another 
chart before that happens. 

Define the chart (u2,v 2 ) (defined for R < R + ) by setting k = B. From (|3"7)l and 
we find that 

/ R \ -jwj ( R 

U2V2 = v-r-) Vr-- 1 



x exp 



R 2 J / 2R 

arctan f 



W Wsjkfi - a 2 



and this provides a regular covering for any subregion with R < R + . 
In the domain of overlap R- < R < R+ the two charts are related by 

\u\f = \u\f , \v\f = \v\Y. (44) 

Figures 2(a), 2(b) are Kruskal-type diagrams which together give a faithful map 
of any subregion covered by a pair of overlapping charts (ui,Vx) and (1/2,^2). The 
maximal analytical e xtens ion of spacetime (|3U)l is analogous to that of a Reissner- 
Nordstrom black hole™ with an exterior event horizon R + and a Cauchy horizon 
R- (cf. Fig. 3). Here the singularity in the interior of the matter distribution is 
barred by the timelike surface R = ^a m where the spacetime extension ends. 

By taking as the maximum value for a white dwarf to be in equilibrium 
(1.4Mq), we have that the domain for the brane tension is given by the lower limit 

771 

In fact, the behavior of the S function does not change when \a\ > lCP 7 ^-. There- 
fore we justify the typical value \a\ = 0.05 adopted in our numerical illustration in 
Sec. 2. 



4. The Hawking Temperature 

In 1975, S. W. Hawking derived - through a semi-classical approach - the thermal 
spectrum of emitted particles by a black hole!^ In this Section we are going to 
follow this same original procedure in order to derive the corrections in the Hawking 
radiation. 
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Fig. 3. Penrose diagram for the spacetime assuming M > M*. The infinite chain of asymptotically 
flat regions I (oo > R > R+) are connected to regions III > R > 7a m i n ) by regions 

II > R > R-). The shaded portion, limited by i? = -fa m i n and R = -f corresponds to the 

interior of the matter distribution. The same dotted line ABCDE illustrates how the matter surface 
evolves in time. Once crossed R— such a surface keeps bouncing between R = ■y and R = fa m i n . 



Let us consider a massless Klein-Gordon field ip in the background defined by 
the exterior spacetime (|30l) . The propagation of such scalar test field is taken to be 
governed by the scalar wave equation 



.g Q0 V a V 6 <p = 0. 

Exploiting the symmetries of the background we seek a solution as 



—Ruii{R*)Y m i( 
Jx 



) exp(— iojt). 



(45) 



(46) 



Substituting this expression for ip, the wave equation is reduced to an ordinary 
differential equation in R* for the modes R u i and given by 



d 2 R u 



dR 



2G N M 3G N M 



f2 



1 - 



R ir\a\R 4 
2G N M 3G N M 2 



R 



AirlalR 4 



)} 



Ruii — 0. 



(47) 
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As r — > ±oo we have that 



d 2 R, 



, u^Rui = => R u i(R*) = exp(±iuR*) 
an* 

and, asymptotically, one can express the Klein Gordon field as 

fi = 4 cxpH^(t - iT )]F m/ (#, 0) (48) 

and 

p 2 = iexpH^+ii*)]^^,^). (49) 
K 

Let us now assume that the source that generates the exterior solution is 
given by a thin shell of a spherically symmetric matter distribution, and the flat 
spacetime inside such distribution is given by 

ds 2 = -dt 2 + dr 2 + r 2 dfl 2 . (50) 

Defining a(t) as the scale factor that the describes the evolution of the matter 
distribution, we impose that the interior metric match the exterior geometry by the 
following equation 

/<fa\2 _ / 2G N M 3G N M 2 \ rdT\* 
~\dt) V a + 47r|cr|a 4 ) \ dt J 

( ( 2G N M , 3GArM 2 \-i/da\2 



a 47r | ct | a 4 

We also define the respective null interior and exterior coordinates by 

V:=t + r, U:=t-r, (52) 

and 

v:=T + R* , u:=T-R*. (53) 

Let us now assume that the null incident rays get into the matter distribution when 
a = ai 3> R+ ~ GAf . Therefore we have that 



and 



On the other hand, 



where 



a/ 47r|cr|a/ z 



1 -> t ~ T. (55) 
«/ ~ £ + iT =*> Vr = vi + n (56) 
k = cti — -R*((ij). (57) 
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When r = 0, we derive the trivial relation between V and U at the center of the 
matter distribution: 

V = t = t/ - (58) 

Let us now consider that the outgoing waves emerge from the matter distribution 
when a = an ~ If to is taken to be the instant in which a = R + , one may 
expand the scale factor aji(t) in Taylor series as 

au(t) ~R + + F(t -t), (59) 

where F is a constant. 

Therefore, from equation (|5ip we have 

T~-Aln(t -t) (60) 

up to first order in (to — f ) , where 



l2G N MR\ - 3GArM 2 /47r | a 

However, from ([57]) we have 



A = 4i? + — 6GnM ■ (61) 



iJ*~Gln(^--l). (62) 

Then we get 

u/7 ~ -(5 In (t - t), 

where 6 = A + G. Therefore, 

C/ // ~ X exp(-^) +^ (63) 

where 

X = -(l + F) , ^ = t -R+. (64) 

However, at the origin of the coordinate system we have Uq = Vq. Therefore the 
relation between the exterior null coordinates is given by 



(-?)• »=-M^) (65) 



v = v + x exp 
where 

vq := ip — k. (66) 
Using (|4"!)|) we expand ipi u i m in terms of tp2uim as 

/>oo 

yiw(m = / [cC' W ;m ex P (-w«/u) - /3u,<uim exp (ia/u)]c&i/, 
Jo 

where a*/ wim and Bu'uim are the Bogolubov coefficients PS Therefore, it is straight- 
forward to show^l that 

\0Lu>ulm\ = irvSlPu'uim]. (67) 
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However, it follows from the orthogonality propriety of ipiuilm an d <p2u;im that 

$3 K'ulmf ~ |/Wm| 2 ] = 1. (68) 

Therefore we obtain that the spectrum of the average number of created particles 
on the u>lm mode is given by 

N u lm = V \(i^ui m \ 2 = jtt^—, — -• (69) 

^— ' exp (zirouj) — 1 

The above result corresponds to a Planckian spectrum with associated temperature 

(r+ - R-) 

t - = W 1 (70) 



where 



R 3 

c = 



3G N M (2R% - M/4tt\ct\)(R + - R_). 

(71) 



{R\ + aR+ + $) L 2R+{2R+ - 3G N M) 

The Hawking temperature depends on the parameters M and |oj. We note that 
in the extremal case we have i? + = i?_ = 3GjvAf/2 implying that — > con- 
tinuously as R+ — > R-. The observation of Hawking radiation could, in principle, 
allows us to test our results for finite \a\. Another feature, which demands a carefully 
analysis, is related to the entropy. We dedicate Sec. 5 to this subject. 

5. Black Hole Thermodynamics and the Quest for a Statistical 
Mechanics Model for the Entropy of Quasi-extremal Black 
Holes 

Motivated by the analysis of energy processes involving black holes BekensteirPSl 
made the remarkable assumption that the entropy of a black hole should be pro- 
portional to the area of its event horizon and formulated a First Law of Black 
Hole Thermodynamics where the surface gravity of the black hole appeared as pro- 
portional (via dimensional fundamental constants) to a temperature. Bekenstein's 
results however did not involve any fund amental principle of statistical mechan- 
ics. Two years later S. W. Hawking J2^J by examining the quantum creation of 
particles near a Schwarzschild black hole, showed that the black hole emits par- 
ticles with a Planckian thermal spectrum of temperature Tbh = k/2tt (in units 
h = c = Gm = Kb = 1) where k = 1/4M is the surface gravity of the black 
hole. This striking result fits exactly in the Bekenstein formula for the First Law of 
black hole thermodynamics, thus validating Bekenstein's proposals and fixing the 
proportionality factor connecting the entropy and the area of the black hole. 

Nevertheless by considering the classical theory of General Relativity, the sin- 
gularity issue still posed an insurmountable barrier to the task of constructing a 
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model for the interior of the black hole and of counting its degrees of freedom, what 
would eventually lead to a definition of entropy. 

The results of the geometrical Black Hole Thermodynamics of Bekenstein are 
recovered in our model for the case of quasi-extremal black holes with an additional 
term connected to the work done by the variation of the brane tension. Let us 
consider a small deviation from the extremal case, with R± = 3G ^ M ± e where e is 
infinitesimal. Neglecting higher order terms in e we have from P(R±) = that 



'3G^M 2 1 



&-kGn\(t 



(72) 



from which the useful relation is derived, 



, e ^(3G|^, M + _i_ (73) 

and which is valid about the extremal configuration. The condition for horizon 
formation is given by e > 0. The equality in this latter relation defines a curve (cf. 




2 

M 

Fig. 4. The parameter space (|<r|, M 2 ): the shaded area corresponds to black hole configurations 
with two horizons; the white area corresponds to configurations with no black hole formation. 
The condition e = defines the limiting curve (|cr*|,Mt) that corresponds to extremal black hole 
configurations. The region that we are considering here is given by sufficiently large values of \<r\, 
and very near the extremal curve on the shaded side. 



Fig. 4) in the parameter space (|cr|,M 2 ) that corresponds to extremal black holes. 
The region above the curve is the region of black holes (with two horizons) while 
the region below defines configurations with no horizon formation and therefore 
no black holes. In this sense, e corresponds to a small deviation from the curve 
(ler* |, Af*) towards the black hole area of the parameter space. In this approximation 
the Hawking temperature is small and given by 
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where we now restore the constants h and Kb- By defining the the outer horizon 
area as A outC r '■= A~kR\ — ^cxtr + 12-7rGjv-M*e^ , we obtain for the quasi-extremal 
case that 

where Eq. ([75)1 has been used. We can therefore associate the horizon area of the 
quasi-extremal black hole with the geometrical entropy 

K B 
AG N h~ 

a result which is in accordance to Bekenstein's definitionPSl Equation (f75j) is an 
extended First Law with an extra work term connected to the variation of the 
brane tension. For deviations with \<j\ = const, we recover the form of the First Law 
for the Schwarzschild black hole in the brane scenario. 

We are now led to tentatively construct a statistical mechanics model for the 
bouncing collapsed matter (in the quasi-extremal case) with an associated partition 
function engendering a thermodynamics which, under certain assumptions, may 
be connected to the geometrical thermodynamics discussed above. The following 
facts about the quasi-extremal configurations are fundamental to our approach. 
The quasi-extremal configurations are basically characterized by the parameter e 
(cf. Eq. (|72p ) that fixes not only the oscillatory motion of the collapsed matter 
but also all the properties of the extended spacetime, in particular determining 
the geometrical thermodynamic variables of the exterior spacetime. Actually e is 
a measure of all the fluctuations about the extremal case, either occurring in the 
dynamics of the oscillating collapsed matter or determining the horizon fluctuations 
of the exterior geometry. 

To proceed let us consider the dynamical equation for the scale factor expressed 
by the constraint (Tl5|) . with the initial conditions d(0) = and a(0) = 1. Once the 
surface of matter distribution is defined by R = ja, we can define the momentum 
(per unity of mass) at the surface distribution as pb. '■= 7 a so that the dynamical 
equation for the scale factor is given by the Hamiltonian constraint 



H 

where 



r Ek + V ( R ) =0 , (77) 



T/ , m 3G N M 2 G N M fc 7 2 
V(R) : =8^W- — + —- 

Expanding ([77)) in a neighborhood of R = = 3Gat-M»/2 (the extremal configu- 
ration) we obtain 



n := ^ 



V\ . 4 S ,„ . /. 40 

2 + 9G%M* 



3G N M, 
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As the brane formulation must approach General Relativity in the low energy limit, 
it is natural to expect the value of |er*| to be sufficiently large. In this instance 
the brane tension satisfies the inequality \Jl/ (SitC^a* |) <C 1, which implies that 
the term proportional to R in ((78)) can be neglected. The constraint (|78l) is then 
approximately given by 

where 



2V2 //, 40 \ 2^/2 40%/2 2 



~ ~ 3G N M* y V 1 9G 2 N Mf 2 ) ~ 3G N M, 27G%Mf ' (80) 

As we are assuming a non-interacting fluid, the interior particles of the matter 
distribution must also oscillate with a frequency u>. We should note that the first 
term in the second equality of (|80|) corresponds to the oscillation of the matter 
distribution in the extremal case, which is thermodynamically a configuration of 
zero temperature, while the second term corresponds to oscillations generated by 
small deviations from the extremal configuration which depend on the parameter 
e. This same parameter is responsible for horizon fluctuations about the extremal 
case which give rise to the geometrical Hawking temperature associated with the 
exterior spacetime, as seen from (|74|) . In analogy with the result (|74|) of the Hawking 
temperature, we make the provisional assumption that the second term in (|80|) arises 
from fluctuations which we denote thermal, connected to a temperature T$ given 
by 

^-■( gadasB? )' 1 <81) 

where a is an adimensional constant. In this sense, from the point of view of statis- 
tical mechanics, the extremal configuration would obviously have a zero partition 
function, being a zero temperature configuration. However this is not the case in the 
quasi-extremal cases for which a partition function can be constructed by quantizing 
the thermal fluctuations appearing (|80|) as we proceed to show. 

Let us define as the number of Planck masses contained in the matter of the 
extremal case, iV* = M*/mp\. The approximated motion of our system can then 
be interpreted as the 1-dim motion of iV* noninteracting oscillators with frequency 
w, the energy levels of which - under a quantization procedure - will be given by 
E n = (n + l/2)7kj. The fluctuations about the extremal configuration present in 
u) and now parametrized with the temperature Tg will engender quantum thermal 
fluctuations that will have a fundamental contribution in the partition function. 

The canonical partition function of the system may then be expressed as 



00 



n=0 



- (n+-)/3fiw]} 



Z = \ > exp ' 



exp(/3u>h) — 1- 
exp(-N*Puh/2), (82) 
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where f3 = 1/KbTs and the third equality results from T$ being small for the 
quasi-extremal case. 

The free energy is given by 

F = -^W^, (83) 

where 1Z is an appropriate constant. By definition the entropy of the system can 
be calculated as a function of the free energy through the relation S = —dF/dTs, 
resulting in 

40y/2N*Kh _ 5V2nN,TZ 
bstat ~ 27K B G%M^ s ~ 3G N M*a C ' ( j 

where use was made of Eq. (|8ip . Therefore, using the relations (|72 p -(|73 p . we obtain 

dS stat ~±(dM+^-d\a\), (85) 

where we have fixed K := K B N, with N = ~ 1.2727V*. We see that N arises 

naturally as the analog of an Avogadro number for the internal matter distribution 
of the extremal case. Also Eq. (|85l) . which is a first law of thermodynamics for the 
bouncing collapsed matter, validates our definition (|81[) of Ts as a temperature. 

By comparing Eqs. (jT5"j) and (|85p we are led to identify Tg with the Hawking 
temperature Tr which was derived in Section 4 for the exterior spacetime, therefore 
fixing a — 1 in (I5TT) . We should remark that the statistical entropy derived in (|84p 
differs from Bekenstein's entropy S geom (1761) by a zero temperature additive constant 
which corresponds to the area of the extremal black hole about which our treatment 
is made, namely 

K B 

+ S stat . (86) 

In fact we can see that the statistical entropy (|84[) arises just from the quantum 
thermal fluctuations with which the partition function is built, in distinction to the 
geometrical entropy of Bekenstein. 

Notwithstanding these striking similarities, a question that can now be posed is 
how the concepts of Hawking temperature and Bekenstein entropy can be connected 
with the above defined entropy and temperature of the interior bouncing matter. 
In fact, in accordance with the geometrical black hole thermodynamics, one might 
argue that the entropy of the black hole is an external variable connected to the 
event horizon boundary of the exterior gravitational field. Therefore the entropy 
of the collapsed matter would thus be irrelevant for physical processes outside the 
black hole. However, as far as quasi-extremal configurations are considered, we can 
suggest a mechanism of how the quantum thermal fluctuations may connect the 
entropy of the collapsed matter to fluctuations of the outer event horizon. To see 
this, let us consider the fundamental frequency Aoj, of a particle of the oscillating 



December 12, 2012 2:7 WSPC/INSTRUCTION FILE text-ijmpd 



20 R. Maier and I. D. Soares 



matter distribution with mass m p i, given by (cf. (|80[1 ) 



Its associated momentum is given by Ap = y / 2hm p iAuj. According to Heisenberg 
uncertainty principle the uncertainty in its localization AR > H/2 and, using (|87[). 
we obtain 



AR > 



1 27G%M*h 



2e V 40 v / 2 



For illustration, by taking M* as the Chandrasekhar mass (A/* = 1.4M ) we find 
Ai? > 0.4127626961 x 10- 14 m 2 /e. On the other hand R + ~ 3G N M*/2 = 3x 10 3 to. 
Therefore, if e < 10 _18 to, namely smaller than 10 17 Zp (where lp is the Planck 
length), we can assure that the scale of the fluctuations Ai? > R + so that the 
quantum thermal fluctuations that give rise to the entropy (|84l) of the collapsed 
matter might be connected to the fluctuations of the event horizon in the quasi- 
extremal case. In this sense, Eqs. ff75|) and (|85]) validate the identification of both 
entropies (up to an additive constant corresponding to the area of the extremal black 
hole) with the proviso that the temperature responsible for the thermal fluctuations 
is identified with the Hawking temperature. 

An extension of our analysis beyond the quasi-extremal approximation will be 
the subject of a future work. In this instance the collapsed matter distribution should 
be strongly correlated and its description as noninteracting oscillating particles 
should be drastically modified. 



6. Conclusions and Final Comments 



In this paper we have examined the gravitational collapse of pure dust in a 
braneworld scenario. By considering a timelike extra dimension of the bulk we show 
that the classical singularity produced by the collapse can be avoided, resulting 
instead in a bouncing nonsingular solution. For the case of the total mass M of 
the collapsing dust greater than a critical value M* (which depends on the brane 
tension) the smooth junction of this solution with the exterior geometry results 
in a spacetime analogous to the Reissner-Nordstrdm geometry, with an exterior 
event horizon and a Cauchy horizon which enclose the collapsed bouncing matter. 
We construct the analytical extension of the spacetime up to the surface of the 
bouncing collapsed matter where the two geometries are smoothly connected. For 
M = M„ we have the extremal black hole configuration with one horizon only. This 
case plays a fundamental role in our discussion of possible entropy definitions for 
the system. 

In the case of M > , the effects of the bulk-brane corrections on the classical 
observational tests of General Relativity were calculated but not included in the 
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paper. However for completeness some of these results are briefly commented in the 
following. To start let us note that the dynamics of General Relativity is recovered 
in the limit |er| — > oo. Concerning the advance of the planetary perihelion, due to 
the brane-bulk corrections in the geometry, the perihelion precession per revolution 
decreases when compared to the predicition of General Relativity, for the present 
case of a timelike extra dimension. On the contrary, for a spacelike extra dim ension, 
the perihelion precession increases. In both cases the observational datcP^ restrict 
the brane tension \a\ to be greater than 10 -8 m~ 2 . 

Analogously for the bending of light rays passing in the neighborhood of a 
spherical massive body, the effects predicted by General Relativity are attenuated 
due to the corrections in the case of a timclike extra dimension. In the case of 
a spacelike extra dimension the deflection a ngle of the asymptotes increases. It is 
worth remarking that the observational data ^ do not impose a new limit for the 
brane tension. In fact, in this case we only need that \a\ > 10~ 12 m -2 for both 
cases of a timelike or spacelike extra dimension. Summarizing the effects of the 
bulk-brane corrections on the classical tests of General Relativity are attenuated 
(increased) when one considers a timelike (spacelike) extra dimension. Although 
the predictions of General Relativity are corroborated by the observational tests in 
the solar system, the corrections of our model may turn out to be significative when 
we consider higher energy and/or curvature scales. 

The central result of our work is contained in Sees. 4 and 5. In Sec. 4 we cal- 
culated (through a semiclassical approach) of the corrected Hawking temperature 
in our model. We also obtain that, analogous to the case of General Relativity, we 
have the prediction of a zero temperature for the case of extremal black holes. The 
calculation of the modified Hawking temperature allowed us to derive, analogously 
to Bckenstein, a geometric entropy that confirms (for quasi-extremal black holes) 
the classical prediction that the entropy is proportional to the area of the event 
horizon, as shown in the first part of Sec. 5. Although the classical black hole ther- 
modynamics introduced by Bekenstein was validated by Hawking's semiclassical 
derivation of the black body thermal emission of a black hole, with a temperature 
proportional the the gravity surface of the event horizon of the black hole, black 
hole thermodynamics always seemed to possess a heuristic character since no basic 
principle of statistical mechanics was used in its derivation. Indeed, the definition 
of the entropy of black holes is still an open problem and we actually refer to it as 
a geometrical entropy. 

In this context, in the remaining of Sec. 5, we have tentatively constructed a 
simple statistical mechanics model for the entropy of the oscillatory collapsed matter 
in the case of quasi-extremal black holes. In fact, by quantizing the internal degrees 
of freedom of the quasi-extremal configuration, we constructed a statistical partition 
function and a free energy - and assuming that the temperature responsible for 
the thermal fluctuations obeys the same form of the Hawking temperature - we 
derive an entropy that satisfies a first law of thermodynamics strikingly similar to 
the geometrical first law for the quasi-extremal black hole. This statistical entropy 



December 12, 2012 2:7 WSPC/INSTRUCTION FILE text-ijmpd 



22 ij. Maier and I. D. Soares 

differs from Bekenstein's entropy by a zero temperature additive constant which 
corresponds to the area of the extremal black hole about which our treatment 
is made. In fact we can see that the statistical entropy arises from the quantum 
thermal fluctuations with which the partition function is built, in distinction to the 
geometrical entropy of Bekenstein. 

We should finally mention that the presence of the Cauchy horizon i?_ in the 
maximal analytical extension of the geometry (|30l) may pose the question of a 
possible instability of the spacetime as measured by geodetic observers.^ We are 
aware that this may constitute a problem to our statistical treatment of the bouncing 
collapsed matter, but not to the geometrical entropy formation of Sec. 6 for the 
exterior spacetime. We intend to approach this issue in a future publication. 
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